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A review of the theoretical activity in Italy in the research field of Hadronic Physics 
is given. Specific focus is put on phenomenological models based on the effective 
degrees of freedom of constituent quarks, on parton distributions in hard processes 
' in the Bjorken limit and on the possibility of linking the two concepts via evolution 

equations. A brief introduction is given also about the socalled generalized parton 
distributions. 



Oh! 



X 



1. Introduction 



The Quantum Chromodynamics (QCD) is the commonly accepted theory 
of strong interactions of elementary particles. It displays two fundamental 
properties: the asymptotic freedom at very short distances and the con- 
I finement of quarks and gluons inside the observed colorless hadrons. While 

the former allows for spectacular precision tests of the theory in the frame- 
work of the socalled perturbative QCD (pQCD), the latter yet cannot be 
explained from first principles. The reason is that the confinement hap- 
pens at the hadronic scale (typically, the nucleon size) where the running 
coupling, as, is so large to forbid any perturbative approach. There are 
several ways to attack the problem: formulating QCD sum rules, solving 
the theory on the lattice, etc.. Here, I will focus on the attempts to explore 
the nonperturbative regime using ansatze induced from the fundamental 
symmetries of QCD itself. More explicitly, this kind of activity, conven- 
tionally included in the socalled Hadronic Physics, can be reformulated as 
the study of confined systems of strongly interacting quarks and gluons. 
From this point of view, the overlap with the know-how of the Nuclear 
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Physics community is very large, since the hadron can be considered as a 
strongly correlated many-body system. 

In the following, I will review the theoretical contributions in Italy to 
this research field. In Sec. 2 I will focus on different formulations of low- 
energy phenomenological models based on the effective degrees of freedom 
of constituent quarks, with emphasys on the covariance and gauge invari- 
ance of such models. In Sec. 3 I will concentrate on the analysis of hard 
processes at a higher scale, in a typical regime of Deep-Inelastic Scatter- 
ing (DIS), with the aim of exploring the (spin) properties of partons inside 
the hadrons and, possibly, of establishing a link between the concept of 
constituent quark (around Aqcd) and current quark (in the DIS regime) 
through the evolution equations. Finally, in Sec. 4, I will consider the so- 
called Generalized Parton Distributions (GPD), a very promising formalism 
which unifies both exclusive and (semi) inclusive hard processes, as it may 
offer a gauge invariant solution to the puzzle of the nucleon spin decompo- 
sition into the spins of its elementary constituents. 

2. Constituent Quark Models 

Constituent Quark Models (CQM) are described by hamiltonians involv- 
ing effective degrees of freedom, the constituent quarks, which are massive 
quasi-particles moving in a background field generated by gluons. Confine- 
ment is realized by a long-range phenomenological potential Vconf, while 
the short-range interaction V"qq is treated perturbatively. The spectrum 
is obtained by diagonalizing the hamiltonian. Since this is not sufficient 
to discriminate among different dynamical models, hadron electromagnetic 
observables are also deduced by computing matrix elements of the electro- 
magnetic current operator between the obtained eigenfunctions. 

A typical example is represented by the socalled hypercentral CQM of 
baryons, where the unperturbed hamiltonian is given by a linear and a 
coulombic contributions in terms of the hyperradius x = yp^ + A^, with 
p, A, the usual Jacobi coordinates of the three-quark system. Such a three- 
body hypercentral potential satisfactorily reproduces the SU(6) gross struc- 
ture of the light-baryon spectrum using just two free parameters. The fur- 
ther hyperfine splittings were interpreted, in a first version of the model, by 
adopting for Kjq the spin-spin part of the hyperfine interaction induced by a 
one-gluon-exchange (OGE) mechanism ^ (employing two new parameters). 
However, despite the introduction of relativistic corrections both of kine- 
matical and dynamical type to the current operator ^, several drawbacks 
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arise. In particular, the wrong ordering in the spectrum of the Roper reso- 
nance = ^ with respect to its positive-parity partner, and the impos- 
sibihty of finding a consistent parametrization for the root-mean-squared 
(rms) charge radius {r'^)ch and the nuclcon (N) electromagnetic form factors 
at the same time. A new version of the hypercentral CQM has been re- 
cently released ^, which implements Voge by adding also an explicit isospin 
dependence (and introducing four new parameters), as suggested also by 
other models of Vqq (sec Sec. 2. 2 in the following). This SU(6)-breaking 
feature is crucial to get a much more accurate description of the spectrum. 

However, the hypercentral CQM is based on a nonrelativistic descrip- 
tion, which is known to be inadequate for constituent quarks that are es- 
sentially relativistic particles. The question arises why these kinds of non- 
relativistic CQM reproduce several hadron properties at the quantitative 
level. A possible answer stems from the socalled General Parametrization 
(GP) method, according to which the (QCD) exact matrix clement of an 
hermitean operator can be represented by a linear parametrization of the 
complete set of spin-flavor operators of three constituent quarks. The key 
observation is that there is a hierarchy in the size of the parameters, the 
nonrelativistic CQM retaining just the first largest ones. Several applica- 
tions of the GP method lead to interesting results The GP method 
considers matrix elements in the hadron rest frame, thus loosing the prop- 
erty of covariance but keeping Lorentz invariance. 

Another series of works consists in covariant approaches based on a 
given hamiltonian with a finite number of degrees of freedom. Along this 
line of investigation, a relativistic quantum-mechanical formulation is only 
possible in terms of unitary representations of the Poincare group. There 
are at least three unitarily equivalent ways to do it: the instant form (IF), 
the light-front form (LF) and the point form (PF) of dynamics. They 
can be distinguished by looking at the different ways of implementing the 
interaction into a free covariant theory while keeping the covariance and the 
Lie algebra of the 10 Poincare generators, the socalled Bakamjian-Thomas 
(BT) construction. The instant form is the usual context in which all 
relativistic field theories, including QCD, are developed. In the following, 
I will review results obtained by the LF and PF descriptions of relativistic 
CQM, respectively. 
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2.1. Front Form Description 

The starting point is the eigenvalue problem for a Poincare- {V-) covariant 
mass operator Mq + V, with Mq the free mass operator. If the interaction 
V is independent of the total momentum Ptot of the system and commutes 
with the LF spin operator, then the cigcnvahie problem can be rewritten 
in terms of an intrinsic relativistic kinetic operator (involving the N con- 
stituent momenta kj such that J^f = 0) ^^'^ of the potential V = TZVTZ^ , 
where TZ is the Mclosh rotation connecting the spins in the (IF) standard 
representation to the LF representation. When calculating electromagnetic 
observables, the current operator J'' must be a conserved current and must 
fulfill the socalled extended V covariance, i.e. it must be P-covariant and it 
must respect parity (P) and time-reversal (T) invariance. Then, the form 
factors can be expressed in terms of matrix elements of only J+ and the 
socalled angular condition is satisfied, i.e. a relation between the current 
matrix elements must hold such that the number of form factors corre- 
sponds to the number of J+ independent matrix elements. If, moreover, 
the Drell-Yan condition for the momentum transfer is satisfied, i.e. = 0, 
the momentum conservation on the light-cone implies, for systems with 
total spin J < ^, a suppression of many-body contributions to the elec- 
tromagnetic vertex, corresponding, in field theory, to Feynman diagrams 
involving nonvalencc qq pairs. Thus, it is legitimate to use the Impulse Ap- 
proximation (lA) for the current operator Jj'^. Introducing parametrized 
form factors for the constituent quarks, a global fit to the spectrum of 
light baryons and to the N and tt elastic form factors is achieved, which 
constrains the parameters in order to make predictions for the transitions 
to the various N* resonances. Using for V the Capstick-Isgur potential 
based on a linear confining term and on an SU(6)-breaking effective gluon- 
exchange contribution, which includes a Coulomb-like part, a spin-spin part 
responsible for the hyperfine splittings, and a tensor part, a successful and 
consistent picture of light hadrons is obtained 

However, the use of Jj'^ breaks the angular condition. This is particu- 
lar evident in the case of the 'N—>- A transition, where three form factors 
(electric, Ge; magnetic, Gm', coulombic, Gc) correspond to four indepen- 
dent matrix elements. Usiially, three out of the latter ones are chosen to 
be independent and the fourth one is forced to fulfill the angular condi- 
tion, but this selection is arbitrary reflecting in very different results for 
the various multipoles contributing to the transition ^. In the case of a 
system with total spin J < 5, there is a perfect correspondence between 
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form factors and Jj^ matrix elements. However, the angular condition is 
still broken, or, equivalently, the rotational covariance of the current op- 
erator is lost. For the example of the N elastic form factors, the G^{Q^) 
can be deduced from different components of J^^ obtaining two different 
results Moreover, it seems not possible to reproduce at the same time 
the two main evidences for the breaking of SU(6) symmetry: the nonvan- 
ishing neutron electric form factor G^{Q'^) and the deviation of the ratio 
Gm(Q^)/G'm((3^) from the value -1.5. Cardarelli and Simula suggest a 
common recipe to these problems. The central idea is that for a system 
with J < 5 in the = Q frame the induced loss of rotational covariance 
implies that the Jf^ matrix elements depend upon the choice of a null 4- 
vector LO perpendicular to the LF null plane .t+ = 0. The constituent quark 
momenta can then be reparametrized as to be all on shell, as demanded by 
the lA. Consequently, the covariant decomposition of the current operator 
contains new spurious structures explicitly depending on lo. Overcoming 
the angular condition problem means selecting appropriate components of 
Jj^ such that the calculated observables do not depend upon spurious con- 
tributions, as they have to. In this way, good results are obtained for the 
electromagnetic observables of pion and nucleon, while the situation is less 
clear in the case of N inelastic transitions because of the scarce amount of 
data available 

For spin-1 systems it is known that also for g"*" = the contribution of 
diagrams involving qq pairs does not vanish, so that the question arises if 
an lA in a LF model with fixed number of degrees of freedom is possible. 
Melikhov and Simula show that this is possible by extending the previous 
analysis to the tensor structure of the deuteron electromagnetic vertex, with 
the interesting outcome that also "good" components for the initial and 
final LF polarization vectors of the spin-1 system must be selected in order 
to get matrix elements free from spurious contributions. 

An alternative solution to the problem of rotational covariance is dis- 
cussed by Lev, Pace and Salme ^. From the fact that LF Lorentz boosts 
are interaction free the wave function is usually split in a cm and an intrin- 
sic part. Similarly, taking matrix elements of the exact current operator 
J'^ on projected N states in the Breit frame with q || z, it is possible to 
construct a current operator, acting only on the intrinsic Hilbert subspace, 
that is extended P-covariant, hermitean, conserved, and fulfills charge nor- 
malization. For all the properties to hold, the choice of the frame is crucial, 
because the rotational covariance around z is sufficient to grant all the other 
features. Since in the LF these rotations are interaction free, the current 
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operator can be built starting from Jj'^, while the angular condition is still 
satisfied. For the example of the deuteron, improvement is obtained for 
both the magnetic and quadrupole momenta, indicating the importance of 
relativistic covariance even at Q"^ = 0; the three electromagnetic form fac- 
tors can be unambiguously calculated, even if the results are still affected 
by model dependence in the choice of the N-N interaction and of the N 
form factors ^. 

As already stressed, the above recipe of building an "intrinsic" covariant 
operator relies on the choice of the Breit frame with q || £, or equivalently 
for qj_ = 0. Therefore, the pair diagrams do contribute to one-loop electro- 
magnetic interactions in a covariant field-theoretical approach. In Ref. 
Simula discusses the issue of rotational covariance by comparing it with 
the previous recipe based on the elimination of spurious contributions. By 
analyzing the form factors of the pion and of the deuteron, he shows that 
the two approaches arc incquivalcnt, since the LF Lorcntz transformation 
from the g"*" = frame to the q_L = one is interaction dependent. The 
claim is that the pair diagrams can be safely neglected in the = frame, 
but they arc important sX, ^ ^ for light hadrons, disappearing in the 
heavy-quark limit. In another paper de Melo, Frederico, Pace and Salme 
analyze the electromagnetic form factor of a pion-like qq bound state by 
building a covariant model of the Bethe-Salpeter amplitude with a symmet- 
ric vertex function in momentum space, deduced from a lagrangian with 
pseudoscalar pion-quark coupling at the one-loop level. By parametrizing 
the Breit frame with an angle a indicating the position of the momentum 
tranfer q in the 'xz plane, they study the relative weight of the spectator 
and pair diagrams with varying a between the boundary positions a = 
(g+ = 0) and a = 90° ((7+ = q), recovering essentially the same results as in 
Ref. However, by rescaling the constituent quark masses to the N masses 
such that the pion-like system becomes a weakly relativistic spin-0 system 
with approximately the deuteron mass, the interesting result emerges that 
the pair diagram can be safely neglected for < 1 and for any a, i.e. in 
any Breit frame with Qy = 0. 

2.2. Point Form Description 

In the PF realization of relativistic quantum mechanics only the four- 
momentum operator P^, generating the space-time translations, contains 
the interaction. The theory is thus manifestly covariant and the current 
matrix elements can be accurately calculated in any frame. It is convenient 
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to introduce the socalled velocity states, defined in such a way that, under 
a Lorentz transformation A, the momenta and the spins of the constituents 
are rotated by the same Wigner rotation Rw; contrary to the LF formal- 
ism, it is then possible to couple orbital angular momenta and spins in the 
same way as it is done in a nonrelativistic theory. By construction, the ve- 
locity states are eigenfunctions of Mq and of the free momentum operator 

~ MqV^. Interactions are introduced in the momentum operator as 

= pi^ + = MV, where the mass operator is M = Mq + V. The 
BT construction is realized if M is a scalar (or, equivalently, V is rotation- 
ally invariant) and commutes with y (i.e., F is a local interaction); the 
Lie algebra of the Poincarc group, then, takes the form of the eigenvalue 
problem M\tp} = m\tp} for a system of mass m made up by N constituents. 

If the current operator satisfies the same Lie algebra, it can be 
shown to be an irreducible tensor operator of the Poincare group. Then, 
a generalized Wigncr-Eckhart theorem implies that if the matrix elements 
are taken on the eigenstates of P'', they can be written as the product of 
covariant objects (the Clebsh-Gordan coefficients of the Poincare group) 
and invariants G^((5^) (the "PF" form factors). In the Breit frame with 
q II z, the Clebsh-Gordan coefficients become the identity, so that the 
G'^{Q^) can be deduced by directly computing the matrix elements of J'' 
in this frame. Recalling the discussion developed in Sec. 2.1, the angular 
condition is here automatically satisfied if is also P— and T— reversal 
invariant; for a system with total angular momentum J there are only 
2J +1 invariant " PF" form factors. 

By projecting the initial and final states in Breit frame upon the corre- 
sponding velocity states with Vi = —Q/2m and Vf = Q/2m, it is possible 
to directly calculate the needed matrix elements of J'* using the elemen- 
tary degrees of freedom. As in the LF case, the lA is adopted but the 
constituents are here considered point-like. The presence of velocity states 
implies that the momentum conservation constraint cannot be written in 
close form. In other words, the impulse q = mvf — mvi delivered to the sys- 
tem is different from the one absorbed by the active constituent. For this 
reason, this lA in PF is called the Point-Form Spectator Approximation 
(PFSA) The role played by the velocity dependent boosts is crucial in 
determining the correct fall-off of the form factors. 

In Refs. 15,16,17 results are shown for the PFSA in a relativistic CQM 
based on the socalled Goldstone-Boson Exchange (GBE) quark-quark in- 
teraction Vgbe, which is dictated by the spontaneous breaking of QCD 
chiral symmetry at low energy. As such, the model is applicable below the 



February 1, 2008 20:1 WSPC/Trim Size: 9in x 6in for Proceedings 



proc 



8 

energy threshold of chiral symmetry restoration, in any case for ^5 3 4 
(GeV/c)^. At variance with Voge, Vgbe exphcitly depends on flavor; its 
peculiar structure is responsible for the very good reproduction of the low- 
lying light and strange baryon spectra using just five parameters. In par- 
ticular, it offers a natural solution to the problem of the correct position of 
the Roper resonance energy level. Covariant calculations of electroweak N 
form factors have been performed in the PFSA obtaining an overall very 
consistent picture ^"^ . The agreement is particularly impressive if it is con- 
sidered that point-like constituent quarks are used and that only the GBE 
model wave functions are input into the calculation, i.e. without introduc- 
ing further adjustable parameters. Comparison with pure nonrelativistic 
calculations indicates that the relativistic boosts play a crucial role even 
at = 0. Moreover, the neutron charge form factor is highly sensitive to 
the details of the GBE model, since it is mostly given by the small SU(6)- 
brcaking components of the GBE eigenfunctions. Only for the magnetic re- 
sponse the results show an evident deviation from data at > 2 (GeV/c)^ 
that leaves room for quantitative improvements beyond the PFSA . 

2.3. Baryon Spectrum and Chiral Symmetry Breaking 

Chiral symmetry breaking is one of the fundamental QCD features that 
must be taken into account when dealing with the low-energy properties of 
hadrons. In this section, the problem is reconsidered replacing a hamilto- 
nian approach with a phenomenological model based on the Nambu Jona- 
Lasinio (NJL) lagrangian. The strong attractive force between quarks in 
the = 0"*" channel induces instabilities in the Fock vacuum of massless 
quarks, breaking the chiral symmetry and generating a dynamical quark 
mass m* related to the vacuum condensate {qq) and of the same order of 
magnitude as the constituent quark mass parameters used in CQM. For in- 
creasing density and temperature the model reproduces the transition from 
the broken phase to the chiral symmetry restoration. In fact, assuming an 
homogeneous, isotropic quark matter inside the baryon, it is possible to 
calculate the radial dependence of the quark dynamical mass m* (r) for the 
flavor a. in a Local-Density Approximation (LDA). It reflects the density 
change when moving away from the baryon center: inside the baryon the 
density (and, consequently, the Fermi momentum) is very large, m*{r) is 
very small and chiral symmetry is restored; getting through the baryon sur- 
face, the density drops while TO*(r) increases until it reaches its constituent 
quark-like asymptotic value, which is determined by the interaction of the 
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quark with the chiraUy broken vacuum. 

The energy of the vacuum is ascribed to the fiUed negative-energy levels 
of the Dirac sea, which put a pressure Pvac on the baryon. Viceversa, inside 
the baryon the pressure j>b, vac is much smaller since the density is very high. 
The equilibrium is given by the counterbalancing pressure p = Pvac — PN,vac 
induced by the quarks inside the baryon. This stability condition fixes the 
size and mass of the baryon According to different choices of the model 
parameters, the range (r^) — 0.76 ^ 0.80 fm is obtained for the N radius 
in very good agreement with the experimental datum Also the baryon 
octet spectrum is fairly well reproduced within a ±5% deviation, even if the 
assumed degeneracy between u and d flavors prevents from distinguishing 
the proton from the neutron and other isospin multiplets. The deviation 
from the baryon decuplet spectrum is, instead, more evident since no spin- 
spin force has been taken into account in the model, while the octet-decuplet 
splitting is usually attributed to the spin-spin interaction. Introducing an 
effective perturbation, based on the spin-spin part of the OGE potential, 
does not substantially improve the situation 

3. Deep-Inelastic Scattering 

Using the optical theorem and the Operator Product Expansion (OPE), the 
socalled Cornwall- Norton moments of the structure functions, that enter the 
hadronic tensor W^" for inclusive DIS on the nucleon, can be expanded as 

M„(g2) = dxx''-''F{x, Q^) = ^ntK(g')l0g(g')] ^^y,2-l ' (1) 

where the Wilson coefficients Cnt are singular in the limit ^ 0, while 
On^ are matrix elements of regular, symmetric, traceless local operators. 
The Cnt can be calculated within pQCD at the given order n; they describe 
the loop corrections that accurately reproduce the observed logarithmic vi- 
olations to scaling in the Bjorken variable x. The On^ contain information 
on the soft physics and are essentially unknown. They are related to non- 
perturbative power corrections through the expansion in twist powers t of 
l/g^, starting from the leading twist t = 2. Both Cnt and On^ implicitly 
depend on a renormalization and factorization scales (which will be under- 
stood in the following), but in any case the OPE gives an automatic and 
rigorous factorization between hard and soft scale physics, even if the proof 
is limited just to inclusive DIS (and e+e~ annihilation). 

The extraction of the socalled higher twists, which in parton language 
represent multiparton correlations, is of delicate importance in determining 
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the full dependence of the hadronic tensor in the preasymptotic region, 
since they affect accurate extraction from data of the leading-twist parton 
distributions, which are interpreted as parton matrix elements of one-body 
operators. In the nc;xt Sec. 3.1 the issue of higher twist extraction will 
be discussed for the unpolarized and polarized structure functions enter- 
ing W^" . In Sec. 3.2 the properties of leading-twist parton distributions 
will be described, with particular emphasys on the socalled transversity 
distribution, which, being a chiral odd object, is unreachable in inclusive 
DIS and represents the missing cornerstone for completing the knowledge 
of the leading-twist spin distribution of quarks inside a hadron. Finally, in 
Sec. 3.3 I will discuss the possibility of linking the parton distributions to 
the wave functions of constituent quarks by properly defining a nonpertur- 
bative, low-scale, phenomenological input to the Altarelli-Parisi evolution 
equations. 

3.1. Structure Functions and Higher Twists 

In Ref. the moments Mn{Q'^) in Eq. (1) are obtained by fitting the 

world data for inclusive electron scattering on protons and dcuterons in the 
DIS kinematics, including data from SLAC and BCDMS, over the range 
1 < < 20 (GeV/c)^ and also for invariant masses lower than 2.5 GeV. 
Since for each Q"^ the whole range in x is needed, it is necessary to use 
phenomenological fits that interpolate the existing data also in the unex- 
plored kinematics and unavoidably introduce a certain amount of uncer- 
tainty. Since the OPE holds for inclusive reactions only, both elastic and 
inelastic contribution have been taken into account, the proton clastic ones 
being analytically calculated in terms of its form factors G^, G^, while the 
convolution formula with a realistic N momentum distribution is adopted 
for the deuteron. In Ref. it is also remarked that moments in the Nacht- 
mann variable ^ are preferable, since they allow to filter out of M„ spurious 
higher-twist contributions coming from power-like target mass corrections 
due to the contributions from O^*? operators with n' ^ n. 

The expansion on the right-hand side of Eq. (1) is approximated up to 
twist t = 6 by using the calculated pQCD result at Next-to-Leading Order 
(NLO) for the leading twist t = 2, and by introducing effective anomalous 
dimensions 7„4, 7„6 for the twists t = 4, 6, respectively, as input parameters. 
The NLO approximation to the leading twist can be further improved in 
two ways. A Sudakov form factor can be included that accounts for the 
Soft Gluon Radiation (SGR) at any order by using resummation techniques. 
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The correction is required by the Wilson coefRcients C„2 being divergent as 
log(n) for X ^ 1, where the cancehation between singularities produced by 
virtual and real gluon loops is lost Alternatively, at the one-loop level 
there are power corrections deriving from vacuum polarization insertions in 
gluon lines, called Infra-Red (IR) renormalons, that cannot be resummed 
in Sudakov form factors; at low they can simulate twist corrections, 
making the extraction of the latter rather ambiguous. 

By fitting the Nachtmann moments on the left-hand side of Eq. (1) with 
the effective anomalous dimension approximation on the right-hand side, it 
is possible to constrain the discussed parameters and extract information 
on the higher twists. The main result is that there is a cancellation be- 
tween the effective twist-4 and twist-6 terms, but each one individually is 
not negligible; this outcome still leaves open the question about a possible 
alternating sign in the twist expansion. It seems, anyway, crucial to simul- 
taneously extract all the twists from data. In fact, a comparison with the 
corresponding analysis obtained by fitting the leading twist with a NLO 
parametrization from Gliick, Reya and Vogt (GRV), indicates that small 
differences at leading twist reflect in larger ones at the much smaller higher 
twists. Most interestingly, the comparison between the calculated leading- 
twist M„ for the deuteron structure functions and the corresponding GRV 
parametrization suggests that the discrepancy of the latter with data can 
be eliminated by enhancing the ratio d{x)/u{x) for the u, d distributions in 
the region x — > 1, to which the higher n moments are mostly sensitive 
This suggestion is consistent with the observed pQCD result | for the ra- 
tio (.t)/_F|' (x) in the elastic limit, against an experimental (questioned) 
datum of J. IR renormalons overestimate the higher- twist contribution 
with respect to the SGR correction, which is preferable also because it is 
analytically calculable. 

The above sketched procedure can be extended to inclusive polarized 
electron scattering to extract Nachtmann moments of the polarized struc- 
ture functions. To this purpose, a new phenomenological parametrization 
of G\{x, Q^) has been developed by interpolating the world data on pro- 
ton targets in the DIS regime for x > 0.02 and 1 < < 50 (GcV/c)^, 
as well as the available results for photo- and electro-production of proton 
resonances. The parameters are constrained to reproduce the experimental 
value of the Gerasimov-DrcU-Hearn (GDH) sum rule. Accordingly, the Q^- 
dependent generalized GDH sum rule is predicted to have a zero-crossing 
point at = 0.16 ± 0.04 (GeV/c)2. For the first moment, the higher 
twists are very small. For higher-order moments, particularly sensitive to 
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the X ^ 1 region, the SGR correction reduces the size of higher twists, 
which, however, remain still relevant for of few (GeV/c)^, at variance 
with the unpolarized case. This suggests that spin-dependent multi-parton 
correlations could have a bigger impact than the spin-independent ones. Fi- 
nally, comparing the trend in the Nachtmann variable ^ of the interpolated 
data for GP(C, Q^) at different low values of and at = 20 (GeV/c)^ 
in the DIS regime, an explicit violation of the Bloom-Gilman duality for 
polarized electroproduction of A on protons is reported. 

3.2. Parton Distribution and Fragmentation Functions 

Assuming factorization between hard and soft physics happening at very 
different scales, the hadronic tensor W^" for inclusive DIS on nucleons can 
be analyzed by the approach based on the diagrammatic expansion. The 
leading-order contribution, both in twist and powers of (aslog((5^)), is given 
by the well-known handbag diagram that involves, in the light-cone gauge 
= 0, the soft quark-quark correlation function 

^a{p) = j dHe-^^ P{P, Si^^iOM^m S) , (2) 

which describes the behaviour of a parton with flavor a and momentum p in 
a target hadron with momentum P and spin S through the hadronic matrix 
element of a leading-twist bilocal operator. Applying the Fierz decomposi- 
tion to $a(p), an expansion in terms of projections $a(a;) onto a complete 
set of Dirac structures F is obtained, where x = p^/P'^ is the light-cone 
momentum fraction of the parton inside the parent hadron for the dom- 
inant light-cone direction "+" induced by the electromagnetic interaction 
(at the considered leading order, the difference between x and Bjorken vari- 
able will be consistently neglected). The T = 7"*", 7"'"75, «a-*+75 projections 
give the leading-twist parton model distribution functions: the momentum 
distribution /i (x) , the helicity distribution gi (x) and the transversity dis- 
tribution hi {x) . Therefore, while three distributions are needed to reach a 
complete knowledge of the leading-twist momentum and spin structure of 
the quarks with respect to a preferred longitudinal direction, only two of 
them, /i and gi, have been extensively studied through measurements of 
the corresponding unpolarized and longitudinally polarized structure func- 
tions, respectively. The reason is that the transversity hi is related to 
helicity flipping mechanisms; hence, it is a chiral-odd object and requires 
the occurrence of another chiral-odd partner in any observable to be mea- 
sured. In inclusive DIS at leading twist, there are no mechanisms providing 
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such partners because chirality is preserved in hard pQCD. Nevertheless, 
hi is a very interesting object. I forward the interested reader to the very 
comprehensive report of Ref. on the hi properties. Here, it is sufficient to 
mention that its integral measures the tensor charge, whose determination 
would give access to chiral-odd operators in QCD and, consequently, to the 
role of chiral symmetry in the hadron structure. Because of the mismatch 
in helicities, hi has a non-singlct evolution disconnected from gluons, hence 
more smooth than the one of gi. In the hadron rest frame, hi and gi would 
be perfectly equivalent because of the rotational symmetry of the system. 
However, the distributions are defined after the system is boosted along a 
preferred longitudinal direction induced by the electromagnetic interaction, 
thus breaking the rotational invariance. Therefore, the difference between 
hi and gi tells us also about the relativistic nature of the dynamics of 
quarks inside the hadron. 

Originally, extraction of hi was proposed by measuring spin asymme- 
tries in Drell-Yan reactions with two colliding transversely polarized pro- 
tons; however, even if the process occurs at leading twist, the probability 
of finding transversely polarized antiquarks is presumable much smaller 
than for quarks. In DIS, semi-inclusive processes need to be considered to 
provide the necessary chiral-odd partner through the quark hadronization 
into one or several detected hadrons. Similarly to Eq. (2), the correspond- 
ing soft correlation function can be parametrized in terms of the socalled 
fragmentation functions, which also have at leading twist a probabilistic 
interpretation and are interesting quantities by themselves since they bear 
witness of how confinement arises. However, the twist analysis reveals 
that for detected spinless hadrons chiral-odd fragmentation functions ap- 
pear only at the twist-three level. The price to pay for having such objects 
at leading twist is to detect baryons, such as the A particle, in socalled 
double-spin asymmetries where both the target and the A are polarized. 
Before commenting on this option, I will discuss the simpler alternative 
of the single-spin asymmetry, where only the target is polarized and final 
mesons (basically pious) are abundantly measured in the detector. The 
key observation is that if memory of the transverse quark dynamics is kept 
inside the various distributions, it is possible to build an azimuthal asym- 
metry related to the socalled Collins angle sin oc x k- St, where Ph is 
the detected hadron momentum with a nonvanishing transverse component 
with respect to the jet axis and k is the fragmenting quark momentum with 
transverse polarization St By properly weighting the single-spin asym- 
metry obtained by flipping the transverse polarization of a nucleon target, 
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it is possible to isolate hi at leading twist by means of the first moment 
of the socalled Collins function H^{z,kT), where z = Pf^ /k~ is the light- 
cone momentum fraction of the fragmenting quark carried by the detected 
hadron. This chiral-odd fragmentation function is given by the probability 
difference of quarks with opposite transverse polarization to fragment into 
an unpolarized hadron; it can be interpreted as a sort of analyzing power 
of the quark polarization into the orbital angular motion of the detected 
hadron. 

The Collins effect has been measured at Fermilab using the reaction 
ppT _> and by the HERMES collaboration using the DIS ep — > e'nX, 
where the longitudinal polarization of the target along the beam direction 
has a small transverse component with respect to the direction of the vir- 
tual photon Because of such an effective small polarization and of an 
asymmetry dilution due to NLO corrections introduced by SGR techniques, 
it came recently as a surprise that the measured asymmetry for all detected 
isospin states of the pion has a remarkable absolute size of as much as 5%. 
At present, there is an ongoing discussion about the issue if factorization is 
under control or if other possible mechanisms could generate a nonvanish- 
ing single-spin asymmetry that could make the hi extraction questionable. 
I refer the interested reader to Ref. and references therein. Here, I will 
just mention the work of Ref. where the Weinberg's idea of a generalized 
time-reversal operator is applied to a flavor multiplet of quarks interacting 
with chiral fields. By mixing up the flavors inside the multiplet on top 
of the usual time-reversal transformation, the internal N chiral dynamics 
generates k^-dependent distributions functions that would be forbidden by 
the standard time-reversal operator, but are perfectly legitimate under the 
new nonstandard one. This would open the possibility of a larger class of 
transverse single-spin asymmetries. Because of the supposed universality 
of such naive T-odd distributions, the same possibility should be open also 
in semi-inclusive lepton DIS, thus giving a possible interpretation of the 
HERMES data. 

A very promising alternative to the Collins effect is represented by the 
single-spin asymmetry generated by the correlation between the transverse 
polarization of the fragmenting quark and the orientation of a tt+tt^ pair 
inside the jet in semi-inclusive production of two leading pions. The effect 
is proportional to the azimuthal angle sin (/) oc P7r+ x • and can be in- 
terpreted as due to the interference between different relative partial waves 
of the pion pair with different phases. By generalizing the soft correlation 
function to this leading-twist analysis reveals that a new class of 
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chiral-odd (and naive T-odd) fragmentation funetions arises, the socalled 
Interference Fragmentation Functions (IFF), that depend in a comphcated 
way upon the kinematics and the "geometry" of the reaction However, 
in DIS a single-spin asymmetry can isolate hi at leading twist through 
H^, which is an IFF sensitive to the transverse component of the pair 
relative momentum The asymmetry is based on a coUinear factoriza- 
tion theorem, where the dependence can be integrated away leading to 
a cancellation of soft gluon radiation effect, at variance with the Collins 
function; an analysis beyond tree level seems, therefore, straightforward. 
The could be extracted, in principle, from the corresponding e+e~ an- 
nihilation. At present, its knowledge relies upon models. In the context of 
the spectator model approximation, the single-spin asymmetry for the DIS 
ep^ e'{TT~^Tr~)X has been computed by modelling the IFF as originating 
from the interference between an s-wave direct tttt production and a p-wave 
decay of the p resonance, leading to an asymmetry as large as 1.5% 

I will now go back to the case of semi-inclusive production of the A 
particle. For the semi-inclusive production of transversely polarized A^^ , a 
double-spin asymmetry depends at leading twist upon the product hiHi, 
where Hi is a chiral-odd (naive T-even) fragmentation function related 
to the probability of a transversely polarized quark to fragment into a 
transversely polarized baryon. A sizeable asymmetry has been observed at 
Fermilab, but the low values of A^ transverse momentum explored prevent 
from using any factorization in the analysis of the cross section, thus posing 
questions on the hi extraction from these data. Moreover, modelling Hi 
is a challenge since the mechanism for the polarization transfer from the 
quark to the A is still unknown. The data collected at Fermilab for the 
unpolarized A^ production {pp — > X) also show a very large asymmetry 
around 20%, while pQCD predicts an almost vanishing effect if the process 
is assumed collincar. Since the A polarization seems roughly independent 
from the target, a possible interpretation is given by a new nonperturbative 
spin effect that produces an asymmetry proportional to the azimuthal an- 
gle sin oc k X • Sa . It is parametrized by the polarized Fragmentation 
Function (polFF) -Dj'^, which is related to the probability of an unpolar- 
ized quark to fragment into a transversely polarized hadron The Dij, 
is naive T-odd but chiral even; therefore, it can be easily detected also in 
A productions with charged-current processes. Assuming a factorization 
theorem and charge conjugation invariance to hold, a fit to the data for the 
pBe A^X and pBe ^ A' X can be achieved with a simple parametriza- 
tion for Dij,. The known trends of a large negative A polarization, increas- 
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ing with Xp = PA/Phcam, of its puzzling flat dependence on A transverse 
momenta beyond 1 GeV/c, and of an almost vanishing A polarization, are 
reasonably reproduced However, the Fermilab data again prevent any 
interpretation in terms of factorization. Therefore, the extracted Df^^ can- 
not be used in the ep e'A^X reaction, where only speculations based 
on different scenarios and different choices of parameters are possible, at 
present 

When the A and A arc longitudinally polarized, new possibilities are 
available for deepening the knowledge of parton distributions Working 
at LO in pQCD, the cross section for the process IN — > I'AX has eight 
helicity components which can be reduced to just four by P-invariancc and 
helicity conservation. Assuming charge-conjugation and isospin invariance, 
linear combinations of these four components for p and n targets allow to 
extract information on the polarized strange quark density gf (x) + g{ {x) 
and on the validity of the hypothesis fl{x) = f{{x) and = g{{x), 

that are not reachable in unpolarized DIS By generalizing the analysis 
also to charge-current processes and e+e~ annihilations leading to A and 
A, various paramctrizations can be considered for the unpolarized Df and 
longitudinally polarized G" fragmentation functions for quark flavor a to flt 
recent data from various laboratories in order to test various scenarios 
according to different assumptions in the flavor SU(3) decomposition. 

3.3. Prom Constituent Quarks to Parton Distributions 

It has already been remarked in Sec. 3.1 that the ratio F^(a;)/F2 (^) 

the limit .x ^ 1 is an observable very sensitive to the symmetry properties 
of the underlying theory, which can influence the parton distributions, in 
particular the ratio d{x)/u(x). When studying the quark distributions, 
usually, the second moment of the valence quark distribution is evolved 
back from a low-energy paramctrization to a low-energy scale /ip where 
the short-range perturbative contribution (related to the gluon and sea 
quarks) is negligible and the valence distribution becomes dominant. Then, 
LO or NLO evolution is applied up to the selected experimental scale and 
the structure functions are deduced in the partonic description of quark 
distributions. Since the parton model emerges in the Bjorken limit because 
of light-cone dominance, it can be developed in the hadron rest frame using 
the light-cone formalism. The low-energy paramctrization of the polarized 
valence quark distribution g™'(a;, /Xq) can then be replaced by the light-cone 
momentum density n\ of the parton with helicity A. The latter can be 
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calculated in the light-front formulation of a rclativistic interacting three- 
body sistem (sec Sec. 2.1). Since, moreover, the usual scale /io turns out 
to be very close to the constituent quark mass, the partons in the hadron 
rest frame can be identified with the constituent quarks. Thus, the low- 
energy input to the evolution equations can be related to the CQM N wave 
functions, that are deduced by satisfying covariance and the Pauli principle. 
In this way, it is possible to test at the high DIS scale the mechanisms 
proposed for SU(6) breaking at low scale. The issue has been debated since 
long time, but only recently in Ref. '^^ a correct treatment of the effects 
related to covariance and Pauli principle have been accounted for. For 
the range 0.3 ^ x < 0.7, where the valence quark picture dominates and 
the approach is meaningful, the ratio F.^ (x)/F2'(^) been investigated 
using the two potentials Voge and Vgbe- Both interactions give asymptotic 
values different from the SU(6) result |, but fail to reproduce the high-a; 
trend of the CTEQ5 parametrization fitting the data, even if the ratio shows 
a rather weak sensitivity to QCD radiative effects introduced by evolution. 
The disagreement emerges only if the effects due to Pauli principle are 
properly included 

The same technique can be applied to the study of the orbital angular 
momentum of the hadron constituents The interest in this observable 
is dictated by the common belief that the N total angular momentum of | 
must be shared among its constituents as i = ^AS -I- Ag + Lq + Lg, where 
AE, Ag are the helicities of quarks and gluons, respectively, and Lq, Lg are 
the quark and gluon orbital angular momenta. A gauge invariant defini- 
tion of Lq and Lg involves the socalled Generalized Parton Distributions 
(GPD), that will be addressed in the next Sec. 4. Here, a low-energy model 
calculation of Lq and AS is obtained by calculating matrix elements of the 

— i(kj X Vi) and {nx — n-\) operators on SU(6) wave functions obtained 
in LF description of the hypercentral CQM described in Sec. 2 Melosh 
rotations of the quark spins generate a nonvanishing Lq at the input scale 
Ho even if the N state is in a 5 wave. The model verifies the N angular 
momentum sum rule, which reads iAE(/z§) = 0.228 + 0.272 = 0.5. 
Then, evolution mixes Tip the two contributions, but the effect of Melosh 
rotations and of high momentum components, generated by Vhyp, survives 
at large x. Evolution produces also Lg{Q'^) ^ 0, even if ig(/Uo) = 0. The 
values obtained in the range 1 < < 20 (GeV/c)^ are compatible with 
the ones obtained by QCD sum rules and lattice calculations 
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4. Generalized Parton Distributions 

As already anticipated in the previous subsection, Generalized Parton Dis- 
tributions (GPD) represent a possible tool to address the orbital angular 
momentum of quarks inside hadrons in a gauge-invariant way. Intuitively, 
GPD arc recovered by considering the handbag diagram of inclusive DIS 
and pushing it to the off-diagonal components in the hadron and par- 
ton legs, as it may happen for the leading-order diagram of Deeply Vir- 
tual Compton Scattering (DVCS) or of hard exclusive Meson Production 
(DVMP). For very large of the entering virtual photon, for large cm 
energy s and small momentum transfer t, a factorization theorem can be 
rigorously proven between the hard electromagnetic interaction and the 
soft physics parametrized in terms of a quark-quark correlation function, 
that generalized the one in Eq. (2). The latter can be expanded in the 
relevant Dirac structure displaying, at leading twist, two unpolarized and 
two polarized helicity-conserving and helicity-flipping GPD. Because of the 
factorization theorem, the GPD formalism represents a unifying framework 
for whole classes of inclusive and exclusive processes. In fact, in the forward 
limit the usual parton distributions are recovered. The first moments in x 
produce the electromagnetic form factors Fi{t), F2{t): if the final hadrons 
are N resonances, the transition form factors are also obtained. From GPD 
spatial transverse distributions can be defined that give information on the 
localization of partons inside hadrons. Finally, at f = the m-th mo- 
ments in X are even polynomials of degree m in the skcwedness parameter 

5, which is related to the mismatch in the longitudinal hadron momentum. 
This constraint, called polynomiality property, is a very general one, since it 
follows from Lorentz invariance. When applied to the unpolarized GPD for 
m = 2, a gauge invariant definition of the quark total angular momentum 
is obtained. 

The experimental extraction of GPD is a very challenging task. Re- 
cently, data have been obtained for the beam spin and charge asymmetry, 
and for DVCS, at the DESY and TJNAF laboratories; these measurements, 
however, allow for mapping GPD only in a very limited part of the phase 
space. Therefore, at present their understanding relies mostly on models. 
In the literature, there are basically two different approaches. The first one 
is a phenomenological construction based on a reduction formula where the 
GPD are related to the parton distributions by factorizing the t dependence 
due to the electromagnetic form factors. The latter approach is based on a 
direct calculation with specific dynamical assumptions. In Ref. '^^ , the chi- 
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ral quark-soliton model is considered, based on the instanton realization of 
the nonperturbative QCD vacuum, and a further proof of its internal con- 
sistency is given by analytically checking that the polynomiality property 
is satisfied. As a byproduct, explicit expressions arc obtained for all the 
coefficients of the m-th order polynomial; in particular, the highest-order 
coefficient gives the m-th moment of the socalled D-term, which gives the 
behaviour of the GPD in the nonvalcnce domain — ^ < x < ^. 

By expanding the N state upon valence light-cone wave functions in 
the Fock space for three partons, the off-diagonal soft correlation function 
can be represented as an overlap of the valence light-cone wave functions 
themselves. Very similarly to Sec. 3.3, the GPD can then be represented 
in terms of the off-diagonal light-cone density matrix, which, in turn, can 
be calculated by solving the eigenvalue problem for some quark-quark po- 
tential in a CQM and then by transforming the solutions to the light-cone 
via the LF boosts The correct covariancc of the approach grants the 
good properties of the GPD, namely recovering the parton distributions 
with the correct support in the forward limit and automatically fulffiling 
the particle number and momentum sum rules. However, since the GPD 
are calculated from valence light-cone wave functions, they are defined only 
in the region ^ < x < 1; therefore, their first moment gives the N electro- 
magnetic form factors only in the limit ^ — > 0. In Ref. the leading- twist 
unpolarized GPD are calculated starting from the l^iyp and T^gbe poten- 
tials described in Sec. 2. No particular difference is obtained with the two 
choices, despite the fact the former is SU(6) symmetric and the latter con- 
tains SU(6)-breaking terms. Remarkably, the unpolarized helicity-flipping 
GPD is completely determined by the effect of Melosh rotations in boost- 
ing the CQM wave functions on the light-cone. A strong t and a weak ^ 
dependence are observed. 
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